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Предложен оптимальный алгоритм формирования радиолокационного изображения поверхности Земли при помощи поляри- 
зационных радиосистем с синтезированной апертурой. Проведен сравнительный анализ полученного и базового алгоритмов. 
Анализ показывает, что предложенный алгоритм позволяет обеспечить более высокую степень соответствия радиолокационно- 
го изображения исходной сцене. Получены математические выражения, позволяющие оценивать потенциальные характеристи- 
ки рассматриваемых систем, в частности, разрешающую способность. 


Радиолокационные системы с синтезирован- 
ной апертурой (РСА) являются эффективным 
средством обзора поверхности Земли для форми- 
рования радиолокационных изображений (РЛИ) с 
высоким разрешением, приближающимся к разре- 
шению оптических изображений. В отличие от оп- 
тических изображений, РЛИ могут быть получены 
независимо от времени суток, состояния облачно- 
сти и осадков, что является существенным преиму- 
ществом РСА. Использование поляризационных 
методов в РСА позволяет значительно повысить 
информативность РЛИ и вероятности распознава- 
ния находящихся на них объектов. 

Ведущие специалисты в области РСА [1,2] от- 
мечают ограничения базового алгоритма формиро- 
вания РЛИ и необходимость в разработке более со- 
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вершенных алгоритмов. Для осуществления поля- 
ризационного анализа в РСА используется метод с 
переключением поляризации зондирующего сиг- 
нала, который приводит [3] к возникновению оши- 
бок формирования РЛИ. В [3] предложен более со- 
вершенный метод одновременного оценивания 
всех элементов матрицы рассеяния. Поэтому воз- 
никает задача оптимизации процесса формирова- 
ния РЛИ, которая и решается в данной работе. 

Рассмотрим процесс функционирования РСА в 
наиболее общей постановке, рис. 1. Предполагает- 
ся, что система движется по некоторой произволь- 
ной траектории, размеры зондируемого фрагмента 
поверхности Земли определяет диаграмма напра- 
вленности (ДН) антенны, которая также считается 
произвольной. 


Технические науки 



Рис. 1 . Зондирование поверхности Земли РСА 


Введем декартову систему координат, в которой 
произвольная точка траектории РСА задается век- 
тором Г, а произвольная точка поверхности Земли 
- вектором §. Будем предполагать, что временная 
зависимость Г(/) - взаимнооднозначная. Это по- 
зволяет использовать один параметр траектории 
(время) вместо двух, подразумевая, что второй при 
этом является известным. 

Процесс формирования РЛИ можно рассма- 
тривать как последовательность двух этапов 
(рис. 2): преобразование зондирующего сигнала 
и, (0 в результате рассеяния на поверхности Земли и 
обработка принимаемого сигнала и,(/), то есть вос- 
становление параметров К(§), характеризующих 
рассеивающие свойства поверхности. Далее эти два 
этапа рассматриваются последовательно. 



Рис. 2 . Этапы преобразования сигнала в РСА 


Зондирующий сигнал с произвольной поляри- 
зацией (полностью поляризованный) считаем уз- 
кополосным. В этом случае он представим в виде 
двухкомпонентного вектора Джонса следующим 
образом 

и,(0 = 11, (Оехр {І(о 0 (}, (1) 


где Щ?) - векторная комплексная огибающая сиг- 
нала, со 0 - несущая частота. 


Сигнал (1), будучи излученным антенной 
(предполагается, что антенна излучает одновре- 
менно обе поляризационные компоненты этого 
сигнала), создает электромагнитную волну, удовле- 
творяющую уравнению Гельмгольца [4]. Решение 
уравнения Гельмгольца для напряженности элек- 
трической составляющей поля в произвольной 
точке г свободного пространства имеет вид 


ЕС, г) = аО(г, 1 )\ ], 


і- 


г-НО 


С(|г-Г(0|), (2) 


где О(т.і) - комплексная ДН, пересчитанная к за- 
данной точке пространства, С - функция Грина 
свободного пространства, а - константа. 

Рассеянное поверхностью Земли поле Е (і,г) в 
борновском приближении описывается следую- 
щим выражением [4] : 


ЕДСг) = п|г>(§С)К(8)Ц 


|§-ГС)|+|г-§| 


хС(|8-Г(0|)С(|г-8|)<28, 


( 3 ) 


где К(§) - матричная дифференциальная функция 
рассеяния; интегрирование производится по всей 
рассеивающей поверхности. Таким образом, рассе- 
ивающие свойства поверхности полностью харак- 
теризуются функцией рассеяния. 


Поле (3), принятое той же антенной, формиру- 
ет принимаемый сигнал, который описывается вы- 
ражением 


ЕД/)=а^0 2 ( 8 ,/)К( 8 )Ц 

5 


1-2 


8-*(0 


хС 2 (|§-Г(0|М§, (4) 


где Ь - константа. 

Рассмотрим процесс обработки сигнала. Мате- 
матическая модель принимаемого векторного сиг- 
нала представлена в виде 

2(0 = и ДО + іі(0, (5) 


где и г (0 - детерминированная составляющая этого 
сигнала, комплексная амплитуда которого опреде- 
лена выражением (4), п(0 - векторный шум, ком- 
поненты которого будем считать некоррелирован- 
ными между собой. 

Задача обработки сигнала (5) - восстановление 
функции рассеяния К(§). Функция рассеяния со- 
держится в принимаемом сигнале в виде параме- 
тра, что позволяет интерпретировать эту задачу как 
задачу оценки параметра. Наличие аддитивной шу- 
мовой составляющей в принимаемом сигнале де- 
лает эту задачу статистической. Решение такой за- 
дачи приводит к алгоритму максимального правдо- 
подобия [3]. Оценка максимального правдоподо- 
бия определяется из условия 

К(§) = шах гѵ(2(0/К(§)), (6) 


где гѵ(2(/)/К(§)) - функционал правдоподобия 
комплексной огибающей принимаемого векторно- 
го сигнала. 


Функционал правдоподобия имеет следующий 
вид 


\ѵ(Ъ(і)) = А ехр 


К 


\іт 


-Ъ(1)\ + [Ъ(І)-Ъ(і)\Аі 


где - спектральная плотность белого шума, Т - 
время наблюдения процесса, А - константа; () + - 
знак эрмитова сопряжения. 

Оценка максимального правдоподобия (6) явля- 
ется решением следующей системы уравнений [5] 

1пи<2,(0) = 0, і,],к = 1,2, (7) 

<И(8) 


где 8/5Щ(%) - функциональная производная. 
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Из (5) следует, что математическое ожидание 
наблюдаемого процесса 


т = И,. (0 = аь\ О 2 ( & /)Щв)Ц ( - 2 


х<? 2 (| 8 -4*Ж 

После ряда математических преобразований 
решение системы уравнений (7) принимает вид 


і(§) = ^|к(§')ч / (§,§Ѵ§, 


1(8) = | Я 2 (8,02401^-2^ 

хС 2 ( |§-Г(0К 
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хП 1-2 
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Вид ПВМН зависит от формы зондирующего 
сигнала, траектории РСА и ДН антенны. Выбирая 
должным образом эти параметры, можно добиться 
выполнения условия (11), т. е. обеспечить высокую 
степень соответствия РЛИ реальной поверхности. 
Таким образом, условие (11) определяет задачу вто- 
ричной оптимизации процесса формирования РЛИ. 

Зондирующий сигнал РСА, как правило, явля- 
ется импульсным. Рассмотрим выражения для ал- 
горитма формирования РЛИ и ПВМН в этом слу- 
чае. Представим (9) в виде 

Т 

І( 8 ) = ]■ 2(014,0, 2)Г*, (12) 

0 

где ТГ ге/ 0,8) “ опорный сигнал. 

При импульсном зондирующем сигнале имеем 

4*&е) = Еіи(*-из;,8), (13) 

л-1 

где Т р - период повторения импульсов. 

^,,/0,8) = 

= О 2 ( 2 , Ф р ,[і - 2 ехр(-2і* |2 - 40 1)- (14) 

I с 


волна означает комплексное сопряжение. 



Рис. 3. Матричный согласованный фильтр 


Выражение (8) описывает результат матричной 
пространственно-временной фильтрации, аппа- 
ратная реализация которой имеет вид, показанный 
на рис. 3. На рисунке ПВСФ - пространственно- 
временной согласованный фильтр (согласован с 
первой или второй поляризационной компонентой 
сигнала). Из выражения (8) видно, что результат 
такой согласованной фильтрации представляет со- 
бой функцию рассеяния, сглаженную функцией 
Ѵ В(2,§’), выражение (10), которую назовем про- 
странственно-временной матричной функцией 
неопределенности (ПВМН). Таким образом, 
ПВМН полностью определяет результат оптималь- 
ной обработки (8). Если эффективная ширина мо- 
дулей элементов этой функции мала по сравнению 
с требуемой разрешающей способностью, а ее вне- 
диагональные элементы пренебрежимо малы, то 

І(2)«К(8)- (11) 


Вид выражений (13) и (14) позволяет произве- 
сти следующую замену переменных: 

Г/ = (-пТ р , і,=пТ р . 

После такой замены ур. (13) можно представить 
в виде функции двух переменных 

- ( І2-Ш1П 

= 24(2,014 4-2 ехр(-2й|2-Г(0|).(15) 

ѵ с ) 

Подставляя выражение (15) в (12), получаем алго- 
ритм формирования РЛИ как операцию двумерного 
преобразования сигнала, описываемую выражением 

1(8) = X ° 2 ' (8Л ) ех Р( -2 ік |2 - 40 1) х 

4 

Т / I I ч 

х| 4 -2 |§ ~ Г(0| Ц, . (16) 

о V с 

Структура выражения (16) позволяет разделить 
процесс формирования РЛИ на два этапа, как это 
показано на рис. 4. Первый этап представляет со- 
бой согласованную обработку комплексной оги- 
бающей импульса матричным согласованным 
фильтром (МСФ), второй этап - когерентная обра- 
ботка с компенсацией фазовой модуляции, возни- 
кающей за счет движения РСА, идентичная для 
каждого поляризационного канала, т. е. собственно 
синтезирование апертуры (СА). Подобную струк- 
туру имеет и базовый алгоритм формирования 
РЛИ. Но в отличие от базового, алгоритм (16): 

• предполагает одновременный поляризацион- 
ный анализ [3] ; 

• является верным не только для линейной, но и 
произвольной траектории РСА; 
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Рис. 5. Г оафики модулей элементов ПВМН 
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• учитывает нелинейное запаздывание огибаю- 
щей по апертуре. 



Рис. 4. Структура алгоритма синтезирования при импульс- 
ном сигнале 


Для рассматриваемого случая импульсного зон- 
дирующего сигнала ПВМН имеет следующий вид 

" г 2(г„ (§) -г„ (§'))' 


х Р(8,8') = Х^ 2 (8)А 2 (8')'і < г 


хехр{2г'А:(г„(8)-г„(§'))}, 


(17) 


где 

г п (&) = \&-Щ-пТ р )\, 

Ц> т (т) = с\ѵХ()ѴІ((-т)(іі 

о 


- матрица известных в радиолокации время-ча- 
стотных функций неопределенности [3]. ПВМН 
может быть использована для определения про- 
странственной разрешающей способности РСА в 
соответствии с выражением 

Лі = /|*„(Ьі')Г<& (18) 

где интегрирование производится по всем возмож- 
ным значениям 8- В этом заключается большая 
практическая важность этого математического 
объекта. 


гг (Л _1 и о ех Р{ 17га(2 }’ если 1 - Т и 
У ра\Ч - ] л ^ гг, 

[ 0, если і > Т і ; 

т, {-іпаі 2 }, если Г <71, 

УрПМ)-) п . -г 

0, если I > 71. 

Здесь Ті - длительность импульса. 

Графики модулей элементов ПВМН для таких 
сигналов в координатах дальность - линейная коор- 
дината показаны на рис. 5 для длины синтезирова- 
ния 120 м. Анализ ур. (18) для этих условий показы- 
вает, что разрешающая способность составляет 3 м. 
На рис. 6 показан один элемент ПВМН и ДН антен- 
ны. Видно, что ПВМН имеет периодический харак- 
тер, и если ДН антенны не будет выбрана должным 
образом, на РЛИ будут присутствовать артефакты. 


ДН антенны 


*и 



Рис. 6. Г оафик модуля элемента ПВМН и ДН антенны 


Рассмотрим пример ПВМН (17) для заданного 
зондирующего сигнала и линейной траектории. В 
качестве векторных зондирующих сигналов, удовле- 
творяющих условию (11), предлагается [3] использо- 
вать сигнал, поляризационные составляющие кото- 
рого имеют линейную частотную модуляцию вверх 
или вниз по частоте. Комплексные огибающие этих 
составляющих описываются выражениями 


Таким образом, в работе предложен усовершен- 
ствованный алгоритм формирования радиолока- 
ционного изображения в радиосистемах с синтези- 
рованной апертурой, основанный на физической 
модели рассеивающей поверхности и отвечающий 
статистическому критерию оптимальности. Этот 
алгоритм позволяет формировать радиолокацион- 
ное изображение с учетом эффектов, которые не 
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учитывает базовый алгоритм что, в частности, по- 
зволяет компенсировать искажения, которые дает 
базовый алгоритм при повышении разрешающей 
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